Generation of generalized coherent states with two coupled Bose-Einstein condensates. 
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We present a scheme to prepare generalized coherent states in a system with two species of Bose- 
Einstein condensates. First, within the two-mode approximation, we demonstrate that a Schrodinger 
cat-like can be dynamically generated and, by controlling the Josephson-like coupling strength, the 
number of coherent states in the superposition can be varied. Later, we analyze numerically the 
dynamics of the whole system when interspecies collisions are inhibited. Variables such as fractional 
population, Mandel parameter and variances of annihilation and number operators are used to show 
that the evolved state is entangled and exhibits sub-Poisson statistics. 
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I. INTRODUCTION 

Combined advances in evaporative cooling techniques 
and magneto-optical trapping made it possible to create 
an atomic Bose-Einstein Condensate (BEC) experimen- 
tally, an important achievement of the last decade. Ini- 
tially predicted by Einstein in 1925 [|, it was produced 
in 1995 from a dilute gas of rubidium atoms 2J. Other 
research groups produced condensates using sodium 
lithium |j] and hydrogen |5|. In a second generation of 
experiments, it was shown to be possible to create dou- 
ble condensates. Such a system can be constructed by 
trapping atoms in two different hyperfine sublevels of 
87 Rb 0. 

Measurements of scattering lengths Q and research on 
the dynamics of the relative phase of the condensates B 
and Rabi oscillations of the two BEC populations [lfj| . 
can be carried out by using a laser-induced Raman tran- 
sition in the Rb experimental setup. The experimental 
production of the first BEC and the analogy between the 
behavior of the coherent matter waves and the electro- 
magnetic ones, encouraged the development of Atom Op- 
tics 1 1 1 1 - Nowadays, quantum optics tools are commonly 
used in the study of the BEC properties. Another conse- 
quence of this analogy is the study of problems already 
explored in quantum optics but mapped into the con- 
text of atomic systems. One example is the generation 
of "Schrodinger cat"-like states (SCS) whose creation in 
quantum optics via dynamical procedures involving non- 
linear interaction, was proposed b y Y urke and Stoler ^2] 
and discussed by several authors [Tj, [Tj, [lj, [l(| • These 
schemes involve Kerr-like couplings and, in general, co- 
herent states are used as initial states. 

When two coupled BECs are analyzed by using a 
many-body Hamiltonian within the two-mode approxi- 
mation (TMA), the terms that describe atomic collisions 
are analogous to a Kerr-like interaction 17]. An addi- 
tional Raman transition could be switched on so there is 
a Josephson-like coupling between the two modes. In 
this case, each mode corresponds to one of the BEC 
species and it is necessary to take into account inter- 



and intraspecies scattering processes. TMA was used 
by several authors to explore the possibility of creating 
quantum superposition states in BECs [LsL fl^ | . Cirac et 
al. [l8| calculated the ground state of the TMA Hamil- 
tonian for various choices of coupling parameters. For 
certain sets of parameter values, the ground state is a 
SCS. Gordon and Savage 0], among others, proposed 
the generation of SCS by exploiting the dynamical evo- 
lution of the system, in a similar fashion as has been 
done in the electromagnetic waves 0. Other aspects 
of BECs recently studied are the entanglement d yna mics 
and the generation of entangled states |2Jj, l2ll I22L l23l . 
A dynamical scheme was proposed by Micheli et al. |22( 
in order to generate a many-particle entangled state. In 
their approach, the entangled subsystems correspond to 
the individual atoms in BECs. 

In the present contribution, we propose the generation 
of the Generalized Coherent State (GCS) in a system 
with two coupled BECs. The GCS was introduced by 
Titulaer and Glauber [24| as a generalization of Glauber's 
coherent state, defined as 

|GCS) = exp(-| 7 | 2 /2) V ^Lexp(iv m )\m).(l) 

Here, the phases v m are functions of index m which en- 
sures Poisson excitation statistics. For the periodic case, 
when v m+ i = v m , GCS can be written as a superposition 
of / coherent states with equal mean value of excitation 
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i 

|GCS) periodic = J2 Ck l^ exp ( i2lTk / 1 )) 

k=l 

Within the TMA we determine the conditions at wich two 
BECs, starting as a product of coherent states, first be- 
come entangled and, later, at certain specific times evolve 
to a product of the vacuum state and a GCS. We also 
show that the phases v m of the created GCS are peri- 
odic, and hence it can be rewritten as a superposition 
of I coherent states. The period I is fixed by both, the 
Josephson-like and nonlinear coupling strengths. We also 
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explore numerically the dynamics of the system with the 
same initial state, in which the interspecies collision pro- 
cess is gradually inhibited. In these situations an exact 
GCS in no longer attained, but the evolved state has in- 
teresting properties such as sub-Poisson statistics, at the 
time GCS would have formed. 

The paper is structured as follows. In Sec.^ we review 
the two-mode approximation, defining the parameters of 
interest in our calculation. Section II I II is reserved for 
the analysis of the necessary conditions to obtain a pure 
GCS. Also, we analyze the possibility of controlling the 
number of coherent states in the created superposition, 
by changing the coupling strengths. We estimate the 
evolution time necessary for the formation of the GCS. 
Section Hvl is devoted to the discussion of feasibility and 
sources of decoherence. Section [V] contains a numerical 
calculation of the dynamics of the system when the col- 
lisions between atoms of different species of BECs are 
inhibited. In Section IVII we summarize our results. 



II. THE TWO-MODE MODEL. 

Our system consists of two atomic BECs of different 
atomic species labeled with suffixes a and 6, in a har- 
monic trap characterized by potentials V a ,b (r). Inter- 
action between atoms a and b are well described if we 
assume only two-body collisions. This can be done by 
considering three different scattering processes: a — a, 
b — b and a — b atomic collisions. We are interested in 
the dynamics of this system when Josephson-like cou- 
pling between species a and b of BECs is switched on. 
The second quantized Hamiltonian which describes our 
system is given by [H E |H H3 



H — H„ + Hi, + H„ 



Hr 



(2) 



where 
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(3) 
(4) 
(5) 



and j = a, b. Here, we have omitted spatial dependence 
in quantum field operators, $ a ,b (^ofe)' w hich annihi- 
late (create) atoms at position r. m is the atomic mass 
and V a< b (r) are the harmonic trap potentials and A a ,b 
are the scattering lengths associated with collisions be- 
tween atoms of the same condensate (intraspecies col- 
lisions). Hamiltonian H a b describes the interaction be- 
tween atoms of different species due to two-body colli- 
sions (interspecies collisions). He is the Josephson-like 



coupling between the modes, 8 being the detuning from 
Raman resonance and CI is the Rabi frequency. 

Following a procedure similar to that described in 
Ref. [3, we obtain the TMA Hamiltonian. The field 
operators are written as *q = cf) a (r) a and = 4>b (r) b, 
4> a .b (r) being the real spatial functions associated with 
each mode and a and b the standard bosonic operators. 
Additionally, we consider here 8 = 0, to obtain the total 
Hamiltonian given by 



H 



with 



HLO a a ! a + hU aa a) a) aa + TiWbPb + hU bb b r b r bb 
+2KU ab a)atib - HX (a f 6 + aS f ) , (6) 
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d 3 rcj) a (r) cj) b (r) . 



4>j{r), (7a) 
(7b) 
(7c) 
(7d) 



The TMA Hamiltonian © can be used in the descrip- 
tion of two different experimental situations. The first 
one is the condensation of sodium, where atoms con- 
dense in hyperfine states localized in two different minima 
of the harmonic trap |2^,|2^|. In this case, Josephson- 
like coupling describes tunneling. In some cases, a good 
approximation is obtained by neglecting the interspecies 
collisions. However, it is more general to assume that 

U a b < U aa = Ubb- 

The second situation is connected with the experi- 
ments of the JILA group with condensation of atoms 
on two different hyperfine 87 Rb levels. In this context, 
the Josephson-like coupling is associated with a laser- 
induced Raman transition between the hyperfine lev- 
els. Reported scattering length values follow the relation 
A a : A ab : A b = 1.03 : 1 : 0.97 HQ. From Eqs.JTEJl and 
(fTcjl it is clear that parameters Uy obey the same rela- 
tions, for a fixed spatial mode function (j) a ^ (r). The lat- 
ter is an important condition if we want to use the TMA: 
as we can see from Eqs.Q, the values of the strengths of 
the Hamiltonian depend on the spatial mode functions 
4> a ,b (r). The approximation is valid only if these func- 
tions remain unaltered and the parameters in each term 
of Hamiltonian © can be considered as constants [47| . 
Several authors use A a = Ab = A a b to simplify theoret- 
ical calculations with the Hamiltonian 2jJ, H3 ■ In 
the TMA Hamiltonian JBJ), this situation corresponds to 

U a b — U aa = Ubb- 

In this article, we assume that U aa + Ubb = 2U a b in or- 
der to extend the analytical solution of the Schrodinger 
equation in [3l| and show how the GCS is exactly gener- 
ated. Notice that this assumption applies for both cases: 
equal scattering lengths approximation (U aa — U a b — 
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Ubb) and for the relation between experimental measured 
scattering lengths (U aa : U a b : Ubb = 1-03 : 1 : 0.97). 
Then, using numerical calculations, we explore the situ- 
ation when U a b < U — U aa = Ubb- In this way, we are 
able to study the effect of the interspecies collision term 
on the dynamics and the transition between two differ- 
ent situations which can be related to the experimental 
contexts of rubidium and sodium (U a b ~ 0) condensates. 



III. GENERATION OF GENERALIZED 
COHERENT STATES. 



In this section, we show how the dynamical evolution 
associated with the TMA Hamiltonian can be exploited 
to produce a product of the vacuum state and the GCS. 
We first assume, reasoning by analogy with BECs in op- 
tical lattices [33, that the system could be prepared as 
a product of coherent states 1^(0)) = |a a ) (g) \a h ) where 
ctj are the amplitude of the state thus |<x, | 2 is the atomic 
population on mode-j. It is demonstrated [s3 that the 
manipulation of the Josephson-like coupling, by chang- 
ing the potential depth between the £ local minima of 
the lattice, produces the state n^l ^)- Another rea- 
son is that the coherent state satisfies the conditions for 
full coherence. In experiments, interference patterns be- 
tween two BECs were observed j2^| and collision-rate 
measurements probed the existence of third-order 
correlations. Although similar patterns could be ob- 
tained if BEC state is described either as Fock or coher- 
ent states 0, Ull 1 studies about decoherence process 
due to three-body losses [3(| supports the assumption 
that the state of a BEC is a coherent state with a well- 
defined phase. Also, phase and spatial dynamics were 
explored including the effect of fluctuations by Sinatra 
and Castin [33 and Ref . [27j ■ Results which are in agree- 
ment with the measure of relative phase between coupled 
condensate (9( were obtained by Li et al. [38( considering 
the initial state \a a ) ® \a.b). 

In this work we shall focus on the macroscopic super- 
position state resulting from the evolution of the sys- 
tem itself. Solving the Schrodinger equation associated 
with Hamiltonian ©, as shown in the Appendix lAl the 
evolved state (h — 1) is given by 



and 



!*(*)) 
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[a (tT [P (t)]™ -itU ab (n 2 -m 2 ) 



Xe -2itU ab nm e -iu, t(n+m) 



(8) 



with 



a(t) 
P{t) 



. . sin (Ait) , . . . 

a a cos (Ait) + 1 (Actb - iOia a ) , (9a) 

M 

, . sin (Ait) . . . . 

a.b cos (Ait) + 1 (Aa a + ujiatb) (9b) 

Ai 
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(N) = 


l"a| 2 


+ l«h| 2 , 



(10a) 

(10b) 

(10c) 
(lOd) 



being N the total excitation number of the system, which 
is a constant of motion and Ai the effective Rabi fre- 
quency. 

We see that the state represented by Eq.lJSJ is an en- 
tangled state and there are only two situations where 

(t)) can be written as a direct product: The first is 
when the interaction parameter U a bt is a multiple of 7r 
and the function exp (— 2inmU a bt) in Eq.JSJ is equal to 
unity. Thus, the disentanglement times associated with 
this first condition depend only on the value of the non- 
linear coupling strength U a b- At these times, \tp (t)) can 
be rewritten as a direct product of new coherent states. 
The second case arises at those times such that either 
a (t) or P (t) is zero, and the evolved state can be written 
as a product of the vacuum state and a superposition of 
Fock states. 

The creation of the GCS given by Eq. lllfl is restricted 
to the times associated with the second situation: if, 
for example, at certain evolution time t e the quantity 
a(t e ) = 0, the evolved state can be written as 



l<K*e)> 



|0) <g> e - ^ ^\ 
|0)<8> |GCS). 



-iUahteTl' 



(11) 



Therefore, the GCS is a special superposition of Fock 
states and obeys Poisson statistics. From Eas. (|9a|l and 
(|9bl) we conclude that the condition when either a (t) or 
/3 (t) is zero can be written as 



i tan (Ait e ) 



iwiOj — Ac*i 
Ai 



(12) 



with i ^ j = a or b, depending on which quantity, a(t) 
or (3{t) goes to zero. Next, we analyze two particular 
choices of the interaction parameter, Ait e , leading to the 
GCS: 



1. At times given by 



i i 2 P +1 
Mt P = — - — 71" 



(13) 



where p is a positive or zero integer, the initial 
states satisfy the relation 



Oti = y {±wi + i\i)aj. 



(14) 



Note that in the particular case u> a — ujb and U aa = 
Uab = Ubb IHIlg, we obtain, from Eas. l|10b|) and 
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: and Ai = A, and Ea. lfT^ is 
ia,-. Therefore, the initial mean 



(|10c|) . that oji = 
reduced to a>j = 
number of atoms in each mode, (n a ) and (fib), must 
be equal n ■ ~ 

A0= (p a 



a a\ — \ a b\ ) with the relative phase 
b corresponding to 5. 



2. At times Aii/c = /c7r, we obtain the condition 
dj = 0. Therefore, it is possible to generate the 
GCS if, for instance, initially all the N atoms are 
condensed in the mode-a and the second hyperfine 
level (mode-&) is used as an auxiliary mode. Thus, 
the atomic population leaves mode-a and returns, 
not as a coherent state but as the GCS. For the 



initial state, |*(0)) 
|*(t fc )> = |GCS>„ < 



N) <8> |0>, 
with 



we find that 



|GCS) 



\yNe 



N,0 



-ik ^fi- 
xe *i 



(15) 



This particular case is interesting because it is pos- 
sible to create the GCS without the necessity to 
"imprint" any initial phase relation between the 
coupled BECs. 

It is possible to rewrite the GCS as a superposition of 
coherent states 25] if the phases given by v n = U a bt e n 2 
on Ea. l|ll|l are periodic. In our context, the necessary 
condition to obtain this kind of "Schrodinger cat" -like 
state is U a bt e = r/s, with r and s integers. This im- 
plies that the interspecies collision strength and effective 
Rabi frequency could also be written as a rational frac- 
tion. When this applies, we can use the discrete Fourier 
transform 39] on Ea. (|llf> . It is straightforward to rewrite 
| GCS) as the superposition 



i-i 



\C(t p )) = J2* { , r n S) \P(te)e- 2 ^), 



(16) 



m=0 



where I is the number of coherent states present in the su- 
perposition. This value is defined by the condition below 



/ = 



2s if r and s are odd, 

s if r is even and s odd or vice versa. 



(17) 



The coefficients aL' ' have the form 



1 /_1 

fe=0 



exp 



-in—k 2 
s 



m 
2m— k 



(18) 



We see that the GCS corresponds to a superposition of 
coherent states with the same mean excitation number 
(|/3| 2 ) and relative phases equal to e -2 ^ 12 ?. Note that 
the number of coherent states in the |GCS) depends on 
the ratio of nonlinear U a b to effective Rabi frequency Ai . 



In Figure^we plot the Husimi quasi-distribution func- 
tion at time t e — for uj a = u>b and U aa — Ubb ~ 
10~ 2 U a b- When U a b and Ai are chosen such that r/s = 
2/3, we obtain three distinguishable packets as shown 
in Fig. [Ha). Modifying the U a b/\\ ratio it is possible 
to achieve a superposition of any number of coherent 
states. For instance, if U a b/\\ =8/5 we obtain superpo- 
sitions of five coherent states as shown in Fig. ^b). If 
Uab/^i = 1/2, we obtain eight packages, Fig. dc). Su- 
perposition of nine coherent states, shown in Fig. GJd), 
is obtained when coupling strengths are set such that 
^ = |. The last plot shows that the different gaus- 
sian packets, associated with different coherent states in 
| GCS), start to merge in phase space at high values of 
I. We also note that in Figs.njVd) the deviation from a 
circular pattern of each gaussian packet in the superpo- 
sition arises from interference between the packets, due 
to their proximity. 

Because the effective Rabi frequency depends on the 
values of traps frequencies u)j and collision parameters, 
all the examples above show that the formation of su- 
perpositions of coherent states in this scenario is highly 
sensitive to changes in these quantities and Josephson- 
like coupling. This means that, by controlling the values 
of scattering lengths, the effective harmonic potential, 
and the coupling between two species of condensates, it 
is possible to "build" a superposition of any desired num- 
ber of coherent states, with a defined number of elements, 
mean excitation values and relative phases. Next, we es- 
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FIG. 1: GCS Q-function on plane (Re [7] , Im [7]) at the first 
purification time t e and iV = 25. at, is given by Eg. 1121 with 
aa = VN(l+i)/2. (a)^ = l;. (b)^ = §; (c)^ = |; 



timate the shortest time interval t e required to obtain 
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the product state |0) ® |GCS). The value of t e depends 
inversely on the effective Rabi frequency Ai (t e cc jf-). 
First, we assume Gaussian spatial mode functions fa (r): 



4>i (r) 



/ 1 



\27rrg 



3 -r 2 /4r 2 



where ro = y/h/2mu>, being u> = uj a = u>t,. Then, us- 
ing the typical physical parameters of Rubidium exper- 
iments, oj = 50 s _1 , m = 1.4 x 10~ 25 Kg and the Rabi 
frequency Q ~ 27r ■ 600 s _1 0, we calculate the value of 
Ai from Ea. i|10c|) . Thus, we obtain t e ps 10~ 3 s. It is 
important to note that t e can be set as short as possible 
by varying the Rabi frequency, Q, of Raman transition. 



IV. DISCUSSION OF FEASIBILITY AND 
DECOHERENCE. 

There are several questions about the feasiblity of the 
CGS arising from the results above. The first one is how 
to set the system in a convenient initial state. From all 
the possibilities suggested by Eq. l|12[l . we conclude that 
the most reasonable initial condition is \a a ) <8> |0). This 
state describes a condensate (in a-mode) and an empty 
auxiliary level (fr-mode) described as a vacuum state. In 
this situation, the imprint of a relative phase between a 
pair of coupled BECs is not necessary. 

Second question is the necessity of an efficient atomic 
population transference. If dccohcrcnce affects the pro- 
cess, we cannot guarantee the formation of the state 
|GCS) <g> |0). Following Ruostekoski and Walls Uj the 
effects of decoherence due to noncondensed atoms on 
BECs shows that purity decays fast, being lower than 
0.2 at U aa t w 0.1. Hence, one must have the interaction 
parameter Ait much smaller than the decoherence time 
scale associated with U aa ~ U a b- 

Once the GCS is created and Raman transition is 
switched off, it is necessary to check the effects of both, 
nonlinear interactions and decoherence. Because the 
fraction of noncondensed atoms is small, we can per- 
form a simple calculation assuming that the interaction 
between those atoms and BEC induce phase-damping 
rather than atomic losses. Thus, we shall consider small 
the effect of decoherence due to condensate feeding and 
depleting. The following master equation for the density 



operator of a-mode, p a , applies 



cquatio 

MM 



df>a 

dt 



H a , f>a 



+ K ({hi, p a } — 2fl a pafla) ■ (19) 



with H a = fud a a'a + hU aa o^ <v dd. It is straightforward 
to calculate the solution of Eq. l|19|l . They resemble the 
solutions for the phase-damped oscillator [43|: 



ram fi\ —ituj a (n—m) —itU a a[n(n — l)—m(m — l)] 

Pa \ l ) — e e 

xe -«t(ra-m) 2 p nm r Q y 



(20) 



If k = 0, we are able to study the dynamics associated 
with atomic intraspecies collisions (nonlinear interaction 



term in H a ), assuming that we create successfully the 
GCS described by Eq.jnj. From Eq.J2U, the density 
matrix is given by 



p a (t) = |GCS(i)) (GCS(i)| 



(21) 



with 



|GCS(*)) 

T(t) 
U(t) 



e 2 



vN exp < —i 



UJ 



7T— + (u„ - Uaa) t 

Ai 



7T— h U aa t. 



(22) 



Nonlinear collisions do not affect the character of the 
state and BEC is still in a GCS, with time-dependent 
amplitude T(t) and phase U(t). From the analysis of 
Sec. IIIII we note that superpositions shown in Fig^ can 
be destroyed as time progresses due to the changes on 
function U{t) defined in Eq. l|22|l . The effect of nonlinear 
interaction after the creation of the GCS could be re- 
duced by manipulation of scattering length A a through 
a Feshbach resonance |44|: controlling the scattering 
length, it is possible to change the value of U aa t so U(t) 
varies smoothly with time. 




0.4 0.6 

Time(7/Ly 

FIG. 2: Tr [pa(t)l as function of time. Light gray line: k — 
O.OlUaa, N = 50; Gray line: k = 0.01U aa , N = 100; Black 
line: k = 0.1U aa , N = 50. 

For k ^ 0, we calculate Tr [p 2 (t)] in order to quantify 
the effects of the reservoir of noncondensed atoms on the 
BEC. We plot this quantity in Fig.|2]for different choices 
of k and number of atoms. From this results it is clear 
that, for large values of k the state is no longer a GCS nei- 
ther a pure state. Additionally, the decay rate is sensitive 
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to changes on k and N, decaying faster when both quan- 
tities increase. Hence, the phase damping is a serious 
limitation for manipulation of the GCS. Another source 
of decoherence is the three-body losses. Measurements 
of density-dependent losses demonstrate that three-body 
recombination is the dominant decoherence mechanism, 
which limits the lifetime and size of BECs [3i|. In or- 
der to study the decoherence process due the three-body 
losses, a master equation is derived by M. Jack j3(| , where 
it is shown that a coherent state is a robust state in the 
limit of large- number of atoms. However, it is also shown 
that the superpositions defined in Ea. (|16[) are sensitive 
to the three-body losses. 

Last concern is the effect of temperature and the ratio 
between Josephson and nonlinear couplings. A careful 
study of this effects on dephasing process was performed 
by Pitaevskii and Stringari |45} | . They shown that coher- 
ence is strongly dependent on the ratio between Joseph- 
son coupling and collisions strength, Ujj/X and also with 
temperature, T/X. One must have control over both ra- 
tios to keep them small in order to keep the phase coher- 
ence. 



V. INHIBITION OF INTERSPECIES 
COLLISIONS. 

In this section, we analyze the effect of the interspecies 
collisions on the dynamical evolution (U a b < U = U aa = 
Ubb)- This is done by solving the Schrodinger equation 
numerically by direct diagonalization of the Hamiltonian 
© in a truncated Fock basis {\n a ,rib}}- In order to 
compare the results for U a b < U with those obtained 
when the condition U aa + Ubb = 2U a b is considered, we 
set the initial states as \a a \ = \ab\ with relative phase 
Acf) — 7r/2. We calculate the fraction of the total atom 
population in mode b, (hb)/N = (wb)/N, and its vari- 
ance (\Ahb\ 2 ) = (n 2 ) — (hb) 2 . The "distance" between 
different states in the Fock basis can be analyzed using 
both (|Arib| 2 ) and the variance of operator b, defined as 
(|A6| 2 ) = <6+6> - (&t)(S) This last relation is useful 
to determine whether a given state can be considered as 
an eigenvalue of fib or b. The Mandel parameter 



Q 



\Ah b \ 



rib 



(23) 



and the linear entropy Sb = 1 — Tib \_p 2 W] are used, 
the first to characterize the statistics and the second to 
quantify the purity of the evolved state of mode 6, re- 
spectively. 

It is convenient to recall some well-known values for 
the definitions written above. For a coherent state (|a)), 
associated with b and b' operators, we obtain 



(|An b | 2 ) 
(|A&| 2 ) 



\a\ 
(hb) 
0, 
0. 



(24) 



indicating a Poisson statistics and that \a) is an eigen- 
state of b. For a Fock state, In), we obtain 



(Tib) 

(|An & | 2 ) 
(|AS| 2 ) 
Q 



= n, 

= 0, 

= n 

= -1, 



(25) 



indicating a sub-Poisson statistics and that \n) is an 
eigenstate of the hb operator. In order to compare these 
values with the numerical results, let us calculate the ex- 
pressions above in the case of equal scattering lengths. 
Using the reduced density operator for mode b extracted 
from Eq.(jHJl, it is straightforward to obtain 



n 6 



(|An 6 | 2 ) 
(|Afe| 2 ) 



1/3 wr 



l/3W| 2 { 



_ -2iV[l-coB (2Ut)] 



}• 



Q = 0. 



(26) 



Except for the Mandel parameter, which is time inde- 
pendent, all the functions depend on |/3 (t)\ 2 , which is the 
mean atom population in mode b, written as 



\[3(t)\ 2 = (|a a | 2 + K| 2 )cos 2 A/, 

- 2 ( Q f Q t ~ a l a b sin 2Xt) , 



(27) 



where |« (*) | 2 + |/3(£)| 2 = N. From Eqs.(j2HJ), we can re- 



cover the result obtained from the analysis in Sec. |TTT] 
of the evolved state. The dynamics depends strongly 
on the relative phase and the initial population of both 
condensates. From the behavior of the partial popula- 
tion, (hb), assuming Acf) — and \a a \ = \a b \, we obtain 
|/3 (t) | = \ab\ and there is no transfer of population be- 
tween the condensates. However, if A</> = ir/2 we can see 
that (hb) = \atb\ [1 — sin (2Xt)] and the system undergoes 
Rabi oscillations with period equals to ir/X. 

We also observe that the variance (|A&| 2 ) is zero at 
times corresponding either to ir/U or when |/3(i)| 2 goes 
to zero. Since the reduced linear entropy is also zero 
at these times, as we discussed in Sec lIIII the variance 
(|A6| 2 ) indicates that mode b is in a coherent state. The 
Poisson statistics remains as time passes, independently 
of the entanglement dynamics of both modes. 

In Fig. |3J we analyze the evolution of the atomic frac- 
tion in mode b, (hb) /N, the Mandel parameter Q and the 
linear entropy, 5 b , for decreasing values of interspecies 
collision strength. We also plot the dynamics of each 
variable associated with the condition of equal scattering 
lengths, shown by a solid gray line. The vertical thick 
gray line indicates the time scale for formation of a GCS, 
t e . The first aspect to be noticed in FigEHa) is a shift in 
the effective Rabi frequency of (hb)/N oscillations with 
decreasing U a b- Also, there is an attenuation of Rabi os- 
cillations if we compare both cases U a b = U and U a b = 0, 
shown in the inset. In particular, there are times at which 
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the transfer of population is suppressed and atoms in 
each condensate are trapped. This "self-trapping" phe- 
nomenon was discussed elsewhere ^3 ■ We want to point 
out that the self-trapping can be associated with inhibi- 
tion of interspecies collision and it is found even at slight 
differences between U and U a b- 



<n>/N 



Q 




FIG. 3: Temporal evolution of population of 6-mode, {hb)/N, 
Mandel parameter, Q, and reduced linear entropy, Sb, for 
|a | = |a&|, A(f> = 7r/2 and j- = 2, cv a = u b and U = U aa = 
Ubb- All quantities are dimensionless. U a b = U (solid gray 
line); U ab = 90%i7 (dashed line); U ab = 80%C/ (solid black 
line); U ab = 75%i7 (black dotted line); U ab = 60%U (black 
thick solid line) ; U a b = 25%U (gray dotted line) and U a b = 
(thick light gray line). Vertical thick gray line indicates time 
of GCS formation for equal scattering lengths. 



The dynamics of the Mandel parameter, FigUfb), 
shows that the subsystem state presents sub-Poisson 
statistics at short times, with Q becoming more negative 
as Uab decreases. Super-Poisson statistics are obtained 
at later times. Our results show that an appropriate 
manipulation of U ao leads the initial coherent state to a 
new one, with sub or super-Poisson statistics depending 
on the evolution time. It is interesting to note that, for 
U a b = 0, the time necessary to reach the minimum of 
Mandel parameter Q is almost the same as that of the 
formation of GCS. 

From our results for Sb, Fig.^c), we see that the en- 
tanglement process is now irreversible. A small change 
of U a b produces an increase in the linear entropy and the 
subsystems are unable to recover purity. In view of this 
result, we conclude that slight changes in our conditions 
for generation of the GCS destroys such a state. The 
variance of b operator is shown in Fig. ^ for the same 
values of U a b as in Fig. EI In all cases, there are no times 
at which (|A6| 2 ) = 0, as in the case of equal scattering 
lengths, and mode b never returns to a coherent state. 



-Q 

< 




FIG. 4: Temporal evolution of (|A6| 2 ) for the same values of 
parameters and U a b as in Figure 0. 



At this point, it is worth recalling that in the sodium 
condensate, A is associated with the tunnelling frequency 
between the two minima of potential, which depends on 
the width of the barrier. We can assume the values used 
by Milburn et al. in order to check the time scales in this 
case. With U approximately 53 s _1 and A ~ 0.37 x 10 3 
s _1 , we estimate Tu ~ 6 x 1CP 2 s and t e ~ 2 x 10~ 3 . 
Again, the time scales for the formation of states with 
sub-Poisson statistics is of the order of milliseconds and 
shorter than the time scale associated with the internal 
collisions. 
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VI. SUMMARY 

Using the TMA Hamiltonian for the description of two 
coupled Bose-Einstein condensates, we demonstrate the 
possibility of creating a generalized coherent state in one 
of the condensate modes. The procedure presented here 
implies only dynamical evolution and requires the prepa- 
ration of BECs in coherent states, which must follow the 
condition given by Eq. (|12f> . The time necessary to obtain 
such a state depends only on the effective Rabi frequency 
Ai, which is a function of Hamiltonian parameters. Also, 
it is shown that the ratio between the collision parameter 
U a b and Ai defines the number of coherent states con- 
tributing to the GCS. For U a b < U aa — Ubb, a new kind 
of non-classical statistics state is created. The analysis 
of fractional population, Mandcl parameter and variance 
of the annihilation operator b shows some interesting dy- 
namical effects associated to this state. Such effects are, 
for instance, a shift of the effective Rabi frequency, some 
temporal regimes with sub-Poisson and super-Poisson 
statistics and irreversible entanglement. 

Acknowledgments 

L. S. likes to thank E. I. Duzzioni, F. O. Prado and 
R. M. Angelo for helpful discussions. The authors also 
thanks to the referee for all the valuable critics. This 
work was supported by FAPESP (Fundacao de Amparo a 
pesquisa do Estado de Sao Paulo) under grants 03 /06307- 
9 and 00/f5084-5 and CNPq (Instituto do Milenio de 
Informacao Quantica). 

APPENDIX A: EVOLVED STATE FOR TMA 
HAMILTONIAN WITH U aa + U bb = 2U ab . 

In this appendix, we calculate a general solution of 
Schrodinger equation associated with Hamiltonian ijBJ by 



means of the unitary transformation 

V (7) = e ^( atS ~ aSt ). (AI) 

With this goal, we rewrite the Hamiltonian J5J using the 
number operator, N — n a + and the unbalance pop- 
ulation operator, An = h a — fib- If U aa + Ubb — %U a b = 0, 
we obtain 

H = lu N + wiAn + UabN 2 - A (a r b + <#) , (A2) 

where loq and uj\ are the quantities defined in Eas. (|10a|) 
and l|10b|l . Using the relations 

V^aV = a cos 7/2 + 6 sin 7/2, 

V^tfV = a 1 cos 7/2 + S f sin 7/2, 

V%V = 6 cos 7/2 - a sin 7/2, 

V r tfV = S f cos 7/2 (A3) 
and choosing the unitary transformation parameter 7 = 
arccos (wi/Ai), we obtain the transformed Hamiltonian 

H v = u N + Ua b N 2 + \ 1 An. (A4) 



The effective Rabi frequency Ai = y/ X 2 + lo\ depends 
on the differences between trap frequencies and collision 
strengths Ujj as can be seen from Eas. l|10ll . It is straight- 
forward to find the time propagator operator 

|*(t)) = Ve~ iAvt V^ |tf (0)) . (A5) 

and, considering the initial state |^(0)) = \a a ) ® \atb), we 
finally obtain the evolved state JSJ) with the quantities 
a(t) and f3{t) given by Eq.©. 
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